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1. 


INTRODUCTION 


1 . 1  BACKGROUND 

Moving-base  gravity  gradiometry  began  in  earnest  dur¬ 
ing  the  late  1960s  when  the  Air  Force  Geophysics  Laboratory 
(AFGL) ,  which  was  then  called  Air  Force  Cambridge  Research 
Laboratories,  sought  to  develop  efficient  approaches  to  map 
the  short-wavelength  features  of  the  earth's  gravity  field 
over  large  geographical  areas.  Since  the  state-of-the-art  of 
gradiometry  has  been  reviewed  on  previous  occasions  (Refs,  2, 

3),  only  recent  activity  in  the  field  will  be  reviewed  here. 

With  steady  progress  of  gravity  gradiometer  system  technology 
demonstrated,  the  Defense  Mapping  Agency  (DMA)  decided  in  1981 
to  establish  a  program  to  produce  a  new  moving-base  gradiometer 
survey  system.  Under  this  program,  Bell  Aerospace  Division  of 
Textron,  Inc.  is  developing  a  Gravity  Gradiometer  Survey  System 
(GGSS)  which  is  based  upon  the  successfully  demonstrated  rotating 
accelerometer  gravity  gradiometer  concept  (Ref.  A). 

The  GGSS  program  is  being  pursued  with  the  goal  of 
gradiometrically  surveying  the  gravity  disturbance  vector  over 
wide  areas  to  less  than  1.0  mGal  rms  error  per  component  and 
at  high  density.  In  addition,  the  very  rapid  data  collection 
rate  and  anticipated  low  cost  per  data  unit  of  surveyed  grav¬ 
ity  information  are  sufficiently  encouraging  to  promise  wide¬ 
spread  use  of  airborne  gradiometers  in  the  future.  The  GGSS 
will  be  capable  of  deployment  in  either  an  aircraft  or  land 
vehicle.  Of  these  two  options,  the  recovery  of  the  surface 
gravity  disturbance  vector  from  an  aircraft-based  GGSS  is  more 
difficult  and  is  the  problem  considered  in  this  report. 


The  GGSS  design  consists  of  a  triad  of  gradiometers 
mounted  on  the  azimuth  (inner)  element  of  a  three-axis  iner- 
tially  stabilized  platform.  As  shown  in  Fig.  1.1-1,  the  three 
instruments  are  symmetrically  distributed  about  the  vertical 
axis,  with  each  gradiometer  inclined  at  the  same  "umbrella" 
angle  a.  Figure  1.1-1  also  identifies  the  u,v,w  instrument 
axes  and  local-level  x,y,z  platform  axes.  Each  gradiometer 
measures  two  elements  of  the  gravity  gradient  tensor  in  the 
instrument  u,v,w  coordinate  frame,  i.e.,  the  cross  gradient 
T  and  the  difference  of  the  two  in-line  gradients,  T  -  T 
Mathematical  properties  of  the  gravity  gradient  tensor  are 
reviewed  in  Appendix  A. 


z 


Figure  1.1-1  GGSS  Geometry  and  Coordinate  Systems 

Earlier  works  (Refs.  5,  6  and  7)  have  considered  the 
use  of  gravity  gradiometer  data  for  identifying  geologic  fea¬ 
tures  to  aid  oil  and  mineral  exploration.  This  report  refines 


the  understanding  of  gradiometer  system  error  sources  and  pro¬ 
vides  a  basis  for  both  designing  gravity  field  surveys  and 
anticipating  the  information  content  (and  omission)  of  the 
data  obtained  therefrom. 


1.2  GGSS  PROGRAM  OBJECTIVES 

The  GGSS  is  being  designed  to  provide  the  following 
accuracies  associated  with  recovery  of  the  short-wavelength 
(less  than  500  km)  portion  of  the  gravity  disturbance  vector 
when  operated  as  indicated  in  Fig.  1.2-1  (Ref.  8): 

•  0.9  mGal  rms  error  for  vertical  gravity 
disturbance 

•  0.18  rms  error  for  deflection  of  the 
vertical  (each  axis). 


600m 

ALTITUDE 


Figure  1.2-1  Nominal  GGSS  Airborne  Survey  Parameters 


JL  4' 


The  wavelength  specification  is  imposed  only  for  initial  test¬ 
ing  to  maintain  reasonable  areal  dimensions.  In  larger  surveys 
longer  wavelength  gravity  information  will  be  obtained.  Using 
both  gradiometer  data  and  long-wavelength  gravity  information 
from  other  sources,  the  goal  for  recovery  of  all  wavelengths 
of  the  gravity  disturbance  vector  is: 

•  1.0  mGal  rms  error  for  vertical  gravity 
disturbance 

•  0.2  s^  rms  error  for  deflection  of  the 
vertical  (each  axis). 

These  larger  values  account  for  errors  in  supplemental  gravity 
field  data  (i.e.,  sources  of  long-wavelength  gravity  informa¬ 
tion  such  as  the  one-deg  by  one-deg  mean  anomaly  data  stored 
in  the  DoD  Gravity  Library  at  the  Defense  Mapping  Agency  Aero¬ 
space  Center) . 


GGSS  ERROR  CONTRIBUTORS 


2.1  DESCRIPTION  OF  ERROR  SOURCES 

GGSS-related  errors  and  survey  limitations  represent 
the  two  fundamental  GGSS  error  source  categories.  The  first 
refers  to  inherent  characteristics  of  the  gradiometer  system, 
while  the  second  applies  to  errors  associated  with  the  manner 
in  which  the  survey  is  operated.  A  discussion  of  the  individ¬ 
ual  errors  in  each  category  is  presented  below.  Quantitative 
treatment  follows. 

GGSS-related  errors  -  Contributors  to  GGSS-related 
errors  include  the  following: 

•  Gradiometer  self-generated  measurement 
noise 

•  Environmentally  induced  errors 

•  Navigation,  attitude  and  attitude-rate 
control  uncertainties 

•  Gimbal ,  vehicle  and  nearby-object  self¬ 
gradient  compensation  limitations. 

Survey  limitations  -  Survey  limitations  which  con¬ 
tribute  to  airborne  GGSS  error  include  sampling  effects,  down¬ 
ward  continuation  and  limited  data  extent.  These  are  each 
summarized  below. 


When  inherently  continuous  data  are  sampled,  aliasing 
occurs.  This  effect  represents  an  error  caused  by  the  misrec- 
ognition  of  frequency  components  higher  than  one-half  of  the 


sampling  frequency.  The  bi-directional  survey  traverse  pat¬ 
tern,  consisting  of  orthogonal  tracks,  is  required  to  minimize 
the  effects  of  cross-track  aliasing  since  no  prefiltering  is 
possible  in  this  direction.  A  well-designed  sampling  regime 
and  data  prefilter  are  required  to  minimize  potential  aliasing 
errors  in  the  along- track  direction. 

Downward  continuation  of  measured  gravity  gradient 
data  from  the  survey  aircraft's  altitude  is  required  to  esti¬ 
mate  gravity  disturbance  quantities  at  the  earth's  surface. 

As  a  result,  an  error  due  to  the  loss  of  the  short-wavelength 
portion  of  the  gravity  signal  occurs  because  the  gravity  gra¬ 
dients  associated  with  highly  localized  disturbance  sources 
attenuate  rapidly  with  altitude. 

Limited  data  extent  reflects  the  loss  of  the  long- 
wavelength  portion  of  the  gravity  gradient  signal  because  the 
far  field  is  not  observed.  Thus,  this  error  source  is  strongly 
influenced  by  the  overall  dimensions  of  the  area  which  can  be 
traversed . 

2.2  DETERMINING  THE  CONTRIBUTION  OF  EACH  ERROR  SOURCE 

The  procedure  selected  to  evaluate  the  relative  im¬ 
portance  of  the  GGSS  error  sources  consists  of  comparing  the 
output  error  when  all  error  sources  are  acting  concurrently 
(i.e.,  ordinary  situation)  with  the  hypothetical  case  in  which 
one  error  source  is  totally  eliminated.  By  repeating  this  com¬ 
parison  for  each  error  source,  the  relative  importance  of  all 
the  error  sources  can  be  established.  Although  this  procedure 
provides  an  understanding  of  the  major  error  contributors,  it 
does  not  offer  insight  into  how  small  refinements  in  survey 
design  should  be  made  because  of  the  nonlinear  fashion  in  which 


the  error  sources  act.  Accordingly,  it  is  also  appropriate  to 
compute  the  sensitivity  of  total  survey  error  to  an  incremental 
change  in  each  error  source  while  maintaining  the  optimality 
of  the  filter  (i.e.,  adjusting  the  filter  to  account  for  geom¬ 
etry  redesign).  The  sensitivity  analysis  results  therefore 
identify  specific  factors  which  provide  the  most  control  of 
gravity  disturbance  vector  recovery  accuracy.  Using  the  meth¬ 
odology  outlined  in  Appendix  B,  the  relative  ranking  and  error 
sensitivities  of  the  GGSS  error  sources  are  computed  by  the 
following  five-step  procedure: 


1. 


2. 


A. 


Determine  the  nominal  rms  recovery  error 
with  all  N  error  sources  present;  the 
notation  for  this  error  is  a 


Compute  the  rms  recovery  error  if  the 

i^^  error  source  is  absent;  this  error 
is  denoted  by  o ' 


.  th 


Compute  the  contribution  of  the  i  “  error 
source  acting  alone  if  it  is  independent 
of  the  other  error  sources  and  linearly 
related  to  the  recovery  error.  This 
"conditionally  independent"  contribution 


is  denoted  by  and  given  by 


=  [o' 


(2-1) 


..th 


Scale  the  variance  of  the  i  error  source 
in  Step  3  to  the  sum  of  variances  of  all 
conditionally  independent  error  sources: 

N 


Percent  contribution  =  100  x 


(2-2) 


Compute  the  sensitivity  of  total  recov¬ 
ery  error  to  a  10  percent  reduction  in 
nominal  error  variance  of  each  source 
or  as  appropriate,  a  10  percent  reduc¬ 
tion  in  the  nominal  value  of  each  survey 
parameter.  Applicable  survey  parameters 


are  the  track  spacing,  aircraft  alti¬ 
tude  and  survey  block  size  which  are 
all  identified  in  Fig.  1.2-1.  This 
error  sensitivity  is  expressed  in  terms 
of  a  percent  change  in  total  error 
variance . 

The  above  procedure  reflects  modifications  to  conven¬ 
tional  linear  error  budget  theory  which  were  chosen  to  account 
for  the  nonlinearity  in  the  altitude  error  term.  For  this 
error  source,  Step  2  above  is  accomplished  by  setting  the  sur¬ 
vey  altitude  parameter  h  (which  explicitly  appears  in  the  vector 
transfer  functions  listed  in  Table  B-1  of  Appendix  B)  to  zero. 

The  effect  of  limited  data  on  total  post-survey  resi¬ 
dual  error  is  analyzed  in  a  slightly  different  fashion  from 
Steps  1-5  for  reasons  of  economy.  The  contribution  of  this 
error  source  alone  is  obtained  by  setting  the  effects  due  to 
GGSS-related  errors,  sampling  and  downward  continuation  to 
zero  and  using  the  optimal  error  covariance  equation: 

Cj.  ,  =  C  -  C  C  ,  (2-3) 

6w,6w  w,w  w,z  z,z  z,w  ’ 

where  the  subscripts  6w  and  z  refer  to  the  residual  error  in  the 
gravity  field- related  quantities  w  and  the  observed  gradiometer 
measurements,  respectively.  Further  discussion  of  these  quanti¬ 
ties  (6w,  z  and  w)  is  given  in  Appendix  B. 

The  solution  to  Eq.  2-3  was  implemented  using  the 
GEOFAST  efficient  collocation  algorithm  (Refs.  9,  10  and  11) 
which  exploits  the  special  structure  of  the  gravity  model 
covariance  matrices  and  uses  to  advantage  the  computational 
efficiency  of  fast  Fourier  transforms.  The  Geodetic  Fast  Esti¬ 
mation  (GEOFAST)  algorithm  is  fully  optimal,  resulting  in 


minimum-variance  of  the  estimation  error  when  appropriate  grav¬ 
ity  models  are  used.  The  rms  surv'ey  error  due  to  the  finite 
extent  of  the  survey  was  computed  using  GEOFAST  for  the  case 
in  which  no  error  sources  other  than  finite  data  extent  were 
acting.  The  contribution  of  the  finite  extent  error,  computed 
in  this  way,  was  combined  with  the  contributions  of  the  other 
errors  discussed  previously  using  Step  4  in  the  outlined  proce¬ 
dure  and  root- sum- square  addition. 


3. 


RESULTS 


This  chapter  discusses  the  results  of  evaluating  total 
post- survey  rms  residual  error  caused  by  the  individual  error 
contributors  identified  in  Section  2.1.  The  results  were  ob¬ 
tained  using  the  methodology  described  in  Appendix  B  and  the 
procedures  outlined  in  the  previous  section.  These  results 
are  estimates  of  the  accuracy  improvement  realized  by  either 
totally  eliminating  or  incrementally  reducing  each  error  source, 
as  discussed  earlier. 

Error  contributions  are  evaluated  using  three  gravity 
field  models:  a  "baseline”  model  which  represents  a  "typical" 
area  that  could  be  surveyed,  an  "active"  model  which  repre¬ 
sents  a  "geographically  rough"  area,  and  a  "North  Texas"  model 
which  represents  characteristics  particular  to  the  planned 
GGSS  test  area.  A  good  representation  of  the  variability  of 
the  gravity  field  is  provided  by  these  three  models.  This 
allows  accurate  evaluation  of  the  expected  range  of  variation 
in  GGSS  performance.  The  relationship  between  these  models 
and  actual  physiographic  features  is  discussed  more  fully  in 
Appendix  B. 


b: 


The  total  rms  error  associated  with  recovery  of  the 
surface  gravity  disturbance  vector  from  GGSS  data  collected  at 
altitude  is  given  in  Table  3-1.  The  sensitivity  of  these  errors 
to  the  three  different  characterizations  of  the  gravity  field 
is  indicated.  These  results  are  also  expressed  in  terms  of 
the  percentage  of  allowable  error  based  on  the  GGSS  performance 
specification  defined  earlier.  Note  that  the  results  presented 
in  Table  3-1  indicate  that  the  0.9  mGal  rms  vertical  distur¬ 
bance  accuracy  goal  is  more  challenging  than  the  deflection  of 

3-1 


TAR  I  K  1-1 


SENSITIVITY  OF  GGSS  SHORT -WAVELENGTH  (LESS  THAN  500  km) 
RECOVERY  TO  GRAVITY  FIELD  MODEL 


GRAVITY  FIELD 
MODEL 

RMS  SURVEY  ERROR 

ERROR  AS  A  PERCENT  OF 

RMS  GGSS  SPECIFICATION 

DEFLECTION  OF 
THE  VERTICAL 
(sec) 

VERTICAL 

GRAVITY 

DISTURBANCE 

(nGal) 

DEFLECTION  OF 
THE  VERTICAL 
(%) 

VERTICAL 

GRAVITY 

DISTURBANCE 

(%) 

Baseline 

0.09 

0.7 

50 

78 

North  Texas 

0.11 

0.7 

61 

78 

Active 

0.13 

0.9 

72 

100 

the  vertical  goal  (0.18  sec  rms).  This  situation  is  typical 
for  every  study  case  addressed  in  these  analyses  and  is  conjec¬ 
tured  to  be  true  in  general.  For  this  reason,  subsequent  dis¬ 
cussions  emphasize  the  ability  of  the  GGSS  to  recover  the 
vertical  gravity  disturbance  quantity.  In  addition,  Table  3-1 
shows  that  gravity  disturbances  in  an  area  characterized  by  an 
active  gravity  field  with  unusually  large  high-frequency  con¬ 
tent  are  more  difficult  to  recover  than  a  typical  gravity  field, 
such  as  would  be  encountered  in  the  North  Texas  area.  The 
active  field  corresponds  to  very  rough  terrain.  Thus,  this 
behavior  is  expected  since  rugged  terrain  causes  the  gravity 
disturbance  field  to  possess  much  more  power  (e.g.,  variance) 
in  the  frequency  range  recovered  by  the  gradiometer.  As  a 
result,  more  energy  is  both  aliased  (for  a  given  survey  track 
spacing)  and  lost  to  downward  continuation  (for  a  given  survey 
altitude) . 


GGSS  error  contributor  rankings  and  the  sensitivities  of 
total  GGSS  error  to  each  error  source  are  provided  in  Table  3-2. 
These  results  are  based  on  a  nominal  set  of  survey  operating 
conditions  and  error  model  parameter  values  as  specified  in 
Table  3-3  and  explained  in  Appendix  B.  The  limited  extent  of 
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TABLE  3-2 


GGSS  ERROR  CONTRIBUTOR  RANKINGS  AND  SENSITIVITIES 
FOR  NOMINAL  SURVEY  CONDITIONS 


GRAVITY  FIELD 
MODEL 

ERROR  SOURCE 

PERCENTAGE  OF  ERROR 
VARIANCE  CONTRIBUTED 
BY  EACH  SOURCE 

ERROR  SENSITIVITY 
( PERCENT 
CHANGE  IN 
ERROR  VARIANCE) 

TOTAL 

RMS 

ERROR 

(mGal) 

GGSS-Relaced  Errors 

3A 

6.5 

Survey  Limitations 

Baseline 

•  Sampling  Effects 

12 

9.2 

0.7 

•  Downward  Continuation 

5 

1.7 

•  Limited  Data  Extent 

49 

1.9 

GGSS-Relaced  Errors 

36 

3.3 

Survey  Limitations 

North 

•  Sampling  Effects 

20 

7.4 

0.7 

Texas 

•  Downward  Continuation 

9 

1.5 

•  Limited  Data  Extent 

35 

3.7 

GCSS-Related  Errors 

37 

4.5 

Survey  Limitations 

Active 

a  Sampling  Effects 

30 

16.4 

0.9 

•  Downward  Continuation 

13 

2.8 

•  Limited  Data  Extent 

20 

1.^ 

TABLE  3-3 

NOMINAL  SURVEY  CONDITIONS  AND  ERROR  MODEL  PARAMETER  VALUES 


PARAMETER 

SYMBOL 

VALUE 

UNITS 

Sampling  Interval 

10 

sec 

Track  Spacing 

^2 

5 

km 

Survey  Altitude 

h 

600 

m 

Survey  Speed 

V 

300 

km/hr 

Survey  Block  Size 

L 

500 

km 

Sensor  White  Noise  Floor 

''k 

55 

E^/Hz 

Sensor  Red  Noise  Level 

'■k 

1.7x10"^ 

E^(Hz) 
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1.  A 


the  survey  dominates  all  other  errors  for  the  baseline  gravity 
model  case.  This  is  due  to  the  fact  that  the  nominal  survey 
block  area  simulated  in  this  study  is  500  km.  As  a  result, 
gravity  wavelengths  longer  than  about  one-half  the  dimension 
of  the  survey  region  (i.e.,  about  250  km  in  this  case)  are  not 
recovered  well.  During  actual  survey  operations,  some  modest 
improvement  in  errors  induced  by  limited  survey  extent  is  likely 
since  more  flight  time  may  be  available  (i.e.,  larger  areal 
dimensions  will  be  possible). 

The  nominal  GGSS  design  parameters,  survey  along- track 
sampling  time,  aircraft  altitude,  aircraft  speed  and  naviga¬ 
tion  performance  were  selected  so  that  their  resulting  contri¬ 
butions  to  total  survey  error  would  be  relatively  small  for 
typical  gravity  field  characteristics.  The  5-km  track  spacing 
was  selected  as  a  compromise  between  short-wavelength  gravity 
field  recovery  and  survey  effort  as  measured  in  hours  of  air¬ 
craft  flight  time.  Note  that  in  an  area  involving  rough  ter¬ 
rain  (active  gravity  field),  a  significantly  larger  error 
contribution  results  from  sampling  effects. 

Table  3-2  also  indicates  that,  regardless  of  the  type 
of  gravity  field  which  is  being  surveyed,  the  greatest  benefit 
to  accuracy  improvement  results  from  an  incremental  reduction 
in  the  track  spacing.  This  phenomenon  is  seen  by  the  large 
error  sensitivities  associated  with  sampling  effects  and  is 
not  surprising  since  decreasing  the  track  spacing  provides  the 
most  direct  mechanism  for  improving  estimates  of  the  high- 
frequency  portion  of  the  gravity  field.  The  findings  reported 
in  Table  3-2  indicate  that  the  gradiometer  instrument  accuracy 


is  sufficient  so  that,  with  careful  attention  to  the  surve’ 


density  and  areal  dimensions,  recovery  of  the  gravity  distur¬ 


bance  vector  will  be  achievable  within  specified  accurac 


levels  and  reasonable  time  limits. 


SUMMARY  AND  CONCLUSIONS 


The  GGSS  is  being  developed  by  Bell  Aerospace  Textron 
to  collect  data  for  accurately  mapping  the  short-wavelength 
features  of  the  earth's  gravity  field  over  large  geographical 
areas  in  a  relatively  quick  and  inexpensive  manner.  Detailed 
consideration  of  all  known  error  sources  which  will  affect  the 
GGSS  indicates  that  the  currently  planned  system  and  survey 
design  is  adequate  to  support  gravity  disturbance  vector  re¬ 
covery  to  within  specification. 


Among  the  important  effects  which  contribute  to  survey 
inaccuracy  are  GGSS  hardware  and  system-related  errors,  sampling 
effects,  downward  continuation  and  limited  data  extent.  Estima¬ 
tion  of  the  post-survey  residual  rms  errors  shows  that  errors 
due  to  these  effects  are  somewhat  sensitive  to  the  characteris¬ 
tics  of  the  underlying  gravity  field.  More  high-frequency 
content  associated  with  a  geographical  area  to  be  surveyed 
(i.e.,  situations  involving  rough  terrain)  leads  to  larger 
survey  errors.  Sampling  effects  and  limited  data  extent  are 
generally  the  most  significant  GGSS  error  sources.  Survey 
parameters  can  be  selected  judiciously  so  that  errors  due  to 
these  effects  are  controlled  without  requiring  excessive  sur¬ 
vey  expense.  With  this  careful  approach  to  survey  design, 
components  of  the  gravity  disturbance  vector  can  be  recovered 
with  better  than  one  mGal  rms  accuracy. 


The  gravity  gradient  tensor  is  symmetric  because  of 
the  commutativity  of  mixed  partial  derivatives  (i.e.,  =  T 
In  addition,  since  the  gravity  field  is  conservative,  Laplace 
equation 


T  +  T  +  T  =0 
XX  yy  zz 


(A-4) 


applies  everywhere  outside  the  earth's  surface.  Thus,  com¬ 
plete  determination  of  f  at  any  point  requires  five  indepen¬ 
dent  measurements  which  mandates  the  need  for  three  appropri' 
ately  oriented  gradiometers  in  the  GGSS. 
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APPENDIX  B 

PERFORMANCE  ANALYSIS  METHODOLOGY 


This  appendix  outlines  the  techniques  chosen  to  analyze 
post-survey  residual  gravity  errors.  A  frequency-domain  ap¬ 
proach  is  used  to  evaluate  the  GGSS-related  and  aliasing  meas¬ 
urement  error  sources.  Due  to  the  nonlinear  response  of  survey 
error  to  errors  caused  by  downward  continuation  and  limited 
data  extent,  a  special  error  allocation  treatment  (which  is 
described  in  Chapter  2)  is  required.  Discussion  of  the  two- 
dimensional  frequency-domain  approach  (Ref.  1)  follows. 

The  problem  of  determining  the  total  residual  error 
in  the  estimation  of  the  gravity  field  from  gradiometer  survey 
data  can  be  formulated  in  general  terms  as  follows.  At  every 
point  X  on  a  planar  earth  surface,  the  estimation  errors  6w 
between  the  true  values  of  the  process  w  and  the  best  estimates 
(in  the  mean-square  sense)  w  are  given  as 


6w(x)  =  w(x)  -  w(x) 


(B-1) 


In  the  above  equation  x  =  (X2^,X2)^  denotes  the  position  of  a 
point  on  the  plane  with  x^  and  X2  measured  in  the  east  and 
north  directions,  respectively.  The  components  of  w  are  any 
collection  of  gravity  field-related  quantities. 

The  estimates  w  are  optimally  obtained  from  a  data 
set  of  gradiometer  measurements  z  which  are  linear  combina¬ 
tions  of  gravity  gradient  quantities  y  corrupted  by  additive 
noise  e: 


2(x)  =  y(x)  +  e(x> ;  X  e  M 


(B-2) 


'v-’a' 

r'l 


N'.vSj 


m 


In  addition,  all  measurement  points  (x)  are  assumed  to  form  a 
rectangular  grid  (M)  and  the  measurement  errors  are  modeled  as 
a  stationary  Gaussian  process  independent  of  the  gravity  field. 
(In  practice,  departures  from  the  grid  are  allowed.  Error 
models  discussed  later  in  the  text  account  for  the  inaccura¬ 
cies  induced  by  such  departures.) 

In  general,  a  linear  estimate  of  w  is  given  by 

w(x)  =  K(x)z  (B-3) 

where  all  available  data  on  the  grid  M  are  used  and  the  esti¬ 
mator  K  is  chosen  so  that  the  variance  of  6w(x)  is  minimized. 
This  estimator  is  determined  from; 

•  The  types  of  data  collected  (e.g.,  grav¬ 
ity  gradients)  and  the  relation  of  these 
data  to  the  quantities  being  estimated 

•  The  physical  characteristics  of  the  sur¬ 
vey  (e.g.,  track  spacing,  time  interval 
between  samples,  etc) 

•  A  statistical  model  for  the  gravity  field 

•  Statistical  models  for  the  measurement 
noise  corrupting  the  data. 

The  analysis  of  the  post-survey  residual  errors  6w  is 
performed  in  the  frequency  domain  using  multi-dimensional  Wiener 
filtering.  The  statistics  of  these  errors  at  any  given  point 
depend  upon  the  relative  position  of  the  point  with  respect  to 
the  measurements  and  are  computed  in  the  form  of  power  spectral 
densities  (PSDs).  The  spectral  density  of  6w  is  the  Fourier 
transform  of  its  auto-covariance  and  is  given  by 


with  and  S2  being  spatial  frequencies  measured  in  the  east 
and  north  directions,  respectively.  The  Fourier  transform  of 
the  optimal  estimator  K  can  be  expressed  in  terms  of  the  spec¬ 
tral  densities  of  the  measurement  and  estimated  processes  by 


K(s)  = 


(B-6) 


Since , 


♦z,w  <£>  =  ♦w.z 


(B-7) 


where  #  indicates  complex  conjugate  transpose,  Eq.  B-4  can  be 


rewritten  as 


^ow,6w  ^w,w  ^w,z  ^z,z  ^w,z 


(B-8) 


Taking  into  account  the  independence  of  the  gravity  field  and 
measurement  errors,  it  follows  that 


<t>  =0 

^w,z  WtX 


(B-9) 


In  addition,  Eq.  B-2  can  be  expressed  in  terms  of  spectral 
densities  as 


(b  =  (b  •!■(]) 

^z,z 


(B-IO) 


Therefore,  substituting  Eqs.  B-9  and  B-10  into  B-8  yields 


The  estimated  variables  w  and  measured  quantities  y 
are  related  to  the  anomalous  surface  potential  through  vector 
transfer  functions  G  and  F  as  follows: 


♦w.w  =  ^  ^ 


(B-12) 


(B-13) 


(B-14) 


where  is  the  PSD  of  the  unsurveyed  surface  anomalous 

^o’  ^o 

potential  obtained  from  a  gravity  field  model,  such  as  those 
discussed  later  in  this  appendix.  Appropriate  expressions  for 
the  vector  transfer  functions  G  and  F  are  obtained  from  th? 
entries  given  in  Table  B-1. 

The  spectral  density  of  the  post-survey  residuals  in 
w  is  obtained  from  a  similar  expression,  namely 


'^6w,6w  ■  -  -  ♦dT^.aT 


(B-15) 


O’  o 


where  the  ensemble  spectral  density  of  the  residual  surface 
anomalous  potential  is  given  by 


♦6T„,«T„  =  1/  I* 


(B-16) 


and  where  0  is  the  spectral  density  matrix  of  the  measure- 
ment  errors.  The  vector  transfer  function  G  must  also  take 
into  account  the  transformation  between  the  gradiometer  instru¬ 
ment  frame  and  the  desired  geodetic  east,  north,  vertical 
(x,y,z)  frame.  As  noted  earlier,  each  of  the  two  gradiometer 


TABLE  B-1 

TRANSFER  FUNCTIONS  FROM  ANOMALOUS  SURFACE  POTENTIAL 


QUANTITY 


SYMBOL 


RELATION  TO 
ANOMALOUS 
POTENTIAL 


TRANSFER  FUNCTION  FROM  . 
ANOMALOUS  SURFACE  POTENTIAL^ 


Anomalous  Potential 
at  Height  h 

East  Component  of 
Gravity  Disturbance 
Vector  at  Height  h 

North  Component  of 
Gravity  Disturbance 
Vector  at  Height  h 

Vertical  Component  of 
Gravity  Disturbance 
Vector  at  Height  h 

East-East  Gradient 
at  Height  h 

North-North  Gradient 
at  Height  h 

Vertical-Vertical 
Gradient  at  Height  h 

East-North 
Gradient  at  Height  h 

East-Vertical 
Gradient  at  Height  h 

North-Vertical 
Gradient  at  Height  h 


BT^/Bx 

n 


aTjj/By 


BT. /Bz 
n 


-27ish 


i2nsj  e 


-2nsh 


i2ns^  e 


-2nsh 


-2ns  e 


■2nsh 


T 

XX 

B^Tj^/Bx^ 

.2  2  -2nsh 
-4n  Sj  e 

T 

yy 

a\/ay^ 

,2  2  -2nsh 
-4n  S2  e 

T 

zz 

a^Tjj/Bz^ 

.2  2  -2nsh 
4ns  e 

T 

xy 

B^Tj^/(BxBy) 

.2  -2nsh 

-4n  SjS2  e 

T 

xz 

a^Tjj/(Bx8z) 

..2  -2nsh 

-i4n  SjS  e 

T 

yz 

3\/Oya2) 

. ,  2  -2nsh 

-i4n  82S  e 

*x,y,z  corresponds  to  an  east,  north,  vertical  geodetic  coordinate  frame 
ts  =  (Sj^  +  82^)^'^^;  i  = 


outputs  measures  a  certain  linear  combination  of  the  gradients 
of  the  potential.  This  combination  is  given  by  the  following 
relation: 


VT/6 

-W3 

>/T/6 

VI/3 

■VI/3 


VI/2  <Jb/6 
-VI/6  V6/6 


-1/2  1/3 


-VI/2 

-VI/6 

0 

VI/3 


V6/6 

■VI76 

■VI/3 


(B-17) 


where  represents  the  output  of  the  gradiometer 

(j  =  1.2  and  i  =  1,2,3)  and  T^.  and 

are  elements  of  the  gravity  gradient  tensor  (see  Eq.  A-3)  in 
an  east  (x) ,  north  (y),  up  (z)  reference  frame. 


The  inverse  Fourier  transform  of  (|»g^  yields  the 
covariance  matrix  of  the  residuals:  ~  ~ 


00  00 

5w,6w  J  *^6w,6w  dsj^  ds. 


where 


<X,S>  =  Xj^Sj^  *  X2S2 


(B-18) 


(B-19) 


is  the  scalar  product  of  the  vectors  x  and  s.  The  rms  values 
of  the  residuals  in  the  components  of  w  are  obtained  by  taking 
square  roots  of  the  diagonal  elements  of  the  matrix  resulting 
from  Eq.  B-18,  evaluated  at  x  =  0. 


Equation  B-16  is  the  fundamental  formula  for  evalu¬ 
ating  post-survey  residual  rms  errors.  From  this  equation, 
spectral  densities  and  covariances  of  the  residual  errors  in 


the  estimation  of  the  gravity  field  from  GGSS  measurements  are 
obtained  after  using  Eqs.  B-15  and  B-18. 

The  spectral  density  of  the  measurement  errors,  whose 
inverse  appears  in  Eq.  B-16,  is  given  by 

*^e,e  “  *^2  ’  (B-20) 

where  0^  ^2  PSDs  of  the  GGSS-related  and 

sampling  effects  measurement  error  sources,  respectively.  The 
measurement  error  spectral  density  matrix  0^^  is  diagonal. 
Entries  along  the  diagonal  of  0^^  are  given  by 

0j^(k,k)  =  mj^  +  nj^  (B-21) 

for  -1/(2tj^)  ^  <  1/(2tj^),  k  =  1,2,..., 6  and  with 

=  t2  V  Wj^  (B-22) 

Hk  =  T2  r^/(V  (B-23) 

where  T2  is  the  track-spacing  parameter,  V  is  the  speed  of  the 
aircraft  carrying  the  GGSS,  Wj^  and  are  the  white  noise  floor 
and  red  noise  level  of  the  GGSS  time  spectra  (Ref.  12)  and 
u»2  =  2ns-^.  In  addition,  represents  the  along- track  sam¬ 
pling  interval  given  by 

T,  =  t  V  (B-24) 

1  S 

with  t^  being  the  time  interval  between  samples.  The  nominal 
values  for  these  parameters  are  given  in  Table  3-3.  (The  val¬ 
ues  given  in  Table  3-3  for  w.  and  r.  relate  to  a  double-sided 


PSD  formulation;  the  values  for  a  single-sided  PSD  formulation 
2  -  2  2 

are  17.5  E  (rad/sec)  and  1.0x10  E  /rad/sec,  respectively.) 


For  computational  purposes,  the  aliasing  contribution, 

(>2«  is  conveniently  obtained  directly  from  the  spectral  density 

of  the  unsurveyed  anomalous  surface  potential ,  (>_  ™  ,  using 

^o  ’  ^o 


the  expression 


=  y'  F(s+ej"^A)F^^s+ej‘^A>#T,  T,  (s+ej‘^A),  (B-25) 

^o~  ■  ■  ■ 

T 

where  A  =  (£,m)  is  a  vector  with  integer  components.  In  gen¬ 
eral,  it  is  sufficient  to  only  use  the  terms  corresponding  to 
the  values  of  A  given  by 


A|A  /  0,  I  Is  +  GJ’^AI I  <  2| Is'l I 


(B-26) 


where  s'  is  the  vector  of  Nyquist  frequencies  associated  with 
the  survey  and  I  hi  I  denotes  vector  magnitude.  In  Eq.  B-25 
0  represents  a  rotation  matrix  determined  by  the  orientation, 
a,  of  the  survey  tracks  relative  to  the  east  direction  and  is 
defined  by 


cosa 

-sina 

0  = 

sina 

cosa 

(B-27) 


J  represents  a  spacing  matrix  determined  by  the  separation  be¬ 
tween  survey  tracks  (12)  and  between  samples  along  a  track 
(x,)  and  is  defined  by 


A  stationary  gravity  field  model  can  be  completely 
specified  by  the  spectral  density  of  the  anomalous  surface 
potential,  4>_  ~  .  Spectral  densities,  covariances  and  cross- 

covariances  of  all  other  geodetic  quantities  can  be  derived 

from  ™  using  frequency-domain  techniques  and  appropriate 
^o’  ^o 

transfer  functions  (see  Table  B-1).  Important  features  of  the 
baseline,  active  and  North  Texas  gravity  field  models  used  to 
obtain  the  results  given  in  this  report  are  provided  below. 

The  baseline  and  active  models  each  consists  of  the 
sum  of  two  third-order  Markov  stationary  random  processes 
along  a  single  track.  An  individual  third-order  Markov  model 
(Ref.  13)  has  a  surface  potential  with  a  spectral  density  of 
the  form 

♦  =  lOn  *  (2ns)2]7/2  ^3^29) 

where  1/Pj^  is  the  characteristic  distance,  Oj^  is  the  variance 
of  the  surface  potential,  and  the  magnitude  of  s  is  given  by 

s  =  (s^  +  .  (B-30) 

For  both  the  baseline  and  active  models 

2 

♦t  t  <s>  =  51!  (B-31) 

(k^ 

with  '  as  defined  in  Eq.  B-29.  The  parameters  of  the  base¬ 
line  and  active  models  were  obtained  by  fitting  to  data  from 
the  North  Atlantic  and  Bonin  Trench,  respectively.  These  param 
eters  are  given  in  Table  B-2.  The  models  have  been  related  to 


land  gravity  data  as  well.  On  land,  the  baseline  model  corre¬ 
sponds  to  the  gravity  field  of  moderate  terrain;  the  active 
model  corresponds  to  very  mountainous  terrain. 


TABLE  B-2 

BASELINE  AND  ACTIVE  GRAVITY  MODEL  PARAMETERS 


The  North  Texas  model  (Ref.  14)  is  a  flat-earth  Attenu- 
ated  White  Noise  (AWN)  gravity  model  (Ref.  15)  which  describes 
the  proposed  GGSS  test  area.  A  single  shell  of  an  AWN  model 
contributes  a  surface  potential  with  spectral  density 


♦  (“>(3)  =  83  d2  ,2  /'“'V 


(B-32) 


which  can  be  viewed  as  resulting  from  a  spherical  shell  at 
depth  D^  below  the  surface  of  the  earth  on  which  the  potential 
is  white  and  such  that  the  surface  potential  has  variance  a^. 
The  complete  North  Texas  AWN  model  consists  of  seven  independent 
shells,  such  that 


with  0  as  defined  in  Eq.  B-32.  The  model  was  fitted  to  all 
worldwide  and  local  gravity-field  information  which  is  appli¬ 
cable  to  the  GGSS  test  site.  The  resulting  (flat-earth)  model 
parameters  are  given  in  Table  B-3. 


TABLE  B-3 

NORTH  TEXAS  GRAVITY  MODEL  PARAMETERS 


SHELL 

(k) 

(km) 

‘^k  2 
(m/s)^ 

1 

2.1 

0.023 

2 

5.0 

0.11 

3 

16 

0.72 

4 

52 

5.8 

5 

161 

23 

6 

861 

70 

7 

2150 

330 

Since  the  active  model  contains  a  substantial  amount 
of  high-frequency  energy,  it  is  not  surprising  that  it  applies 
to  areas  with  strong  physiographic  features  such  as  mountain 
ranges,  escarpments  and  ocean  trenches.  Alternatively,  the 
baseline  model,  which  contains  a  moderate  amount  of  both  high- 
frequency  and  low-frequency  energy,  best  characterizes  areas 
containing  more  typical  physiographic  features  such  as  plains, 
shelf s  and  basins.  The  North  Texas  model  is  an  intermediate 
type  model  which  is  appropriate  for  the  physiography  existing 
in  northern  Texas.  For  comparison  purposes,  the  rms  values  of 
the  gravity  anomaly  and  deflections  of  the  vertical  for  each 
of  these  models  are  listed  in  Table  B-A. 


Gravity  Anomaly 

mGal 

- 1 

51.1 

- 1 

112.5 

31.3 

Deflection  of  the  Vertical 
(each  axis) 

1 

sec 

7.6 

I 

16.8 

4.7 

East-East  Gradient 

E 

20.7 

60.6 

13.0 

North-North  Gradient 

E 

20.7 

60.6 

13.0 

Vertical-Vertical  Gradient 

E 

33.8 

99.0 

21.2 

East-North  Gradient 

E 

12.0 

35.0 

East-Vertical  Gradient 

E 

23.9 

70.0 

North-Vertical  Gradient 

E 

23.9 

70.0 
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